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ABSTRACT 

Observation by the CANGAROO-III stereoscopic system of the Imaging Cherenkov Telescope has 
detected extended emission of TeV gamma rays in the vicinity of the pulsar PSR B1706— 44. The 
strength of the signal observed as gamma-ray-like events varies when we apply different ways of 
emulating background events. The reason for such uncertainties is argued in relevance to gamma-rays 
embedded in the "off-source data" , that is, unknown sources and diffuse emission in the Galactic plane, 
namely, the existence of a complex structure of TeV gamma-ray emission around PSR B1706— 44. 
Subject headings: gamma rays: observation — pulsar: individual (PSR B1706-44) — diffuse radiation 



I 1. INTRODUCTION 

PSR B1706— 44 is a young pulsar with a high spin- 
down luminosity, one of the gamma-ray pulsars detected 
by the EGRET instrumen t of the Compton Gamma 
Ray Observatory (Thompso n et al."71 [l992). Detection of 
point-like gamma-ray emission of TeV gamma-rays was 
reported, using data of the CANGAROO-I telescope, at 
the position of the pulsar with a fl ux of ^ fO~^^ cm~ ^ 
s~^ for gamma-ray energy >lTeV ( Kifune et al. Ifl995l ). 
which was not pulsed in modulation with the pulsar 
^ period, and is thus likely, if the detection is valid, to 
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be from a compact pulsar wind nebula associated with 
PSR B1706-44, similar to the case of the Crab neb- 
ula. This result was followed by the Durham Mark 6 
telescope, which detec ted a flux of ~ 4 x 10 ~^^ cm~^ 
s-i at £; > aOOGeV (iChadwick et al. Ifl99l . Obser- 
vations with the CANGAROO-II telescope preliminarily 
reported the detection of a gamma-ray signal, which ap- 
peared to be somewhat broader than the point-spre ad 
function, ~ 0.2° of the telescope (|Kushida et al.ll2C)03[ ). 

In order to reduce background events and to improve 
the accuracy of the arrival direction, stereoscopic obser- 
vations with a system of multiple Imaging Air Cherenkov 
Telescopes (lACTs) have now been put into operation as 
"second generation lACTs", such as H.E.S.S. and VER- 
ITAS, as weU as CANGAROO-III, and in addition, the 
MAGIC telescope of a large aperture of 17m diameter. 
The sensitivity to TeV gamma rays has been dramati- 
cally improved to uncover a number of TeV gamma ray- 
sources, more than 70 at the time of t he 21st Intern a- 
tional Cosmic Ray Conference in 2007 1 (|Hinton 2007f ). 

However, H.E.S.S. reported a null detection of point- 
like emission from the vicinity of PSR B1706— 44, setting 
an upper hmit of 1.4 x 10 ~^^ cm~^ for E > 350GeV 
(|Aharonian et al.~l l2005af ) at the center of the field of 
view. CANGAROO-III failed to confirm the signal, 
and obtained an upper limit of 5 x 10~ ^^ cm~^ s~^ for 
E > 600GeV on the point-like emission (jTanimori et al.l 
I2OOI . 

A number of Galactic sources of TeV gamma rays have 
been discovered through a sc an survey of the Galacti c 
disk by the H.E.S.S. Group (jAharonian et al. 11200 59), 
and it is considered that pulsar wind nebulae (PWN) 
constitute one major class of the sources. A consider- 
able number of those sources are found to exhibit ex- 
tended emission, with its central position displaced from 
the pulsar position in some cases. It may also be the 
case for PSR B1706— 44 that the pulsar is accompanied 
by a PWN of extended TeV gamma-rays, or that some 
unknown object of TcV gamma-ray emission exists in its 
vicinity without having any clear counterpart in other 



2 



R. Enomoto et al. 



wavelengths, such as X-rays. 

We still have a long way to go before the TeV gamma- 
ray signal from lACT will mature to become completely 
free from the background of cosmic ray particles, of which 
intensity outnumbers TeV gamma-rays by orders of mag- 
nitude. In order to reveal gamma-ray events as a promi- 
nent signal, it is required to subtract background events 
from the observed data and, for this purpose, to deduce 
the background events from a dataset which is considered 
not to contain gamma-ray events. 

The H.E.S.S. scan survey of the Galactic disk has 
shown a broad, diffuse emission in the Galactic plane 
near the Galactic Center, when integrated over a solid 
angle of 10"'^ sr, having an intensity as large as the flux 
from the Cr ab pulsar nebula with i ts spectrum as hard 
as ^ E~^-^ (jAharonian et al. ]l2006f ). This fact warns us 
that TeV gamma-ray emission may exist where no promi- 
nent objects are seen, and we have to be careful when we 
use events from an area in the Galactic disk to estimate 
the background. 

Various ways for inferring the background events as 
well as for monitoring the variance of the observation 
conditions are used in different observations. For exam- 
ple, the observations of CANGAROO-I, -II, and -III con- 
sist of two kinds of separate observations of ON-source 
and OFF-source runs. The OFF-source observation, of 
which direction is selected to be more than 30° away 
from the source position, is utilized as the dataset for 
emulating the background in ON-source data. On the 
other hand, the H.E.S.S. group adopts the method of so- 
called "wobble" and/or "ring background", in which the 
datasets for emulating background events are taken from 
"off-source" directions of ^ 0.5° away from the source 
position within the common field of view of the "ON- 
source observation" . 

Presented in this paper is the result of CANGAROO- 
III observations of PSR BI706— 44, obtained through 
analysis by utilizing two different methods for estimat- 
ing the cosmic-ray background by (a) the "OFF-source 
observation" , which is directed away from the Galac- 
tic plane, and (b) the "wobble" and the "ring back- 
ground" method, which uses, for emulating background, 
the events coming from the directions along a ring- 
shaped region within the field of view of the ON-source 
run that contains the observation target near its center. 

2. CANGAROO-III STEREOSCOPIC SYSTEM 

The CANGAROO-III stereoscopic system consists of 
four lACTs located near Woomera, South Australia 
(31°S, 137°E). Each telescope has a 10-m diameter of 
a segmented reflector, consisting of 114 spherical mir- 
rors made of FRF, i. e., fiber-reinforced plastic material 
(|Kawachi et al.ll2001h . each of 80 cm diameter, mounted 
on a parabolic frame with a focal length of 8 m. The total 
area for light collection is 57.3 m^. The first telescope, 
Tl, w hich has been operated as the CANGAROO-II tele- 
scope ()Enomoto et al.ir2002al) . is not presently in use be- 
cause of its smaller field of view than the others and 
a deterioration of the reflectivity of the plastic mirrors. 
The second, third, and fourth telescopes (T2, T3, and 
T4) were used for the observations reported in this pa- 
per. The camera and electronics sys tem for T2, T3, an d 
T4 and other details are given in iKabuki et all (|200l . 
The telescopes are located at the east (Tl), west (T2), 



south (T3) and north (T 4) corners of a diamon d shape 
with its sides of '-lOOm (|Enomoto et al.ll2002bl ). 

3. OBSERVATIONS 

The observations of PSR B 1706— 44 were carried out 
during the period between Jul 11 and 19 in 2004 and 
between Apr 14 and Jun 15 in 2007. During the ON- 
source runs in 2004, the center of the field of view 
of the telescope was set at right ascension a—257.A° 
and declination 5 = —44.5° [J2000], i.e., the posi- 
tion of PSR B170 6-44. In 2007, the "wobble mode" 
(jPaum et al.lll997l ) was adopted, i.e., the pointing posi- 
tion of each telescope was displaced from PSR B1706— 44 
in declination by ±0.5°, and was changed every 20 min- 
utes to be set alternatively at -1-0.5° or at —0.5° from 
S — —44.5°. It is noted that we have a merit, in the 
wobble mode observation, to enlarge the effective field 
of view of the telescope and to average the response of 
the photomultiplier camera which fluctuates from pixel 
to pixel. 

Listed in Table [1] are the directions pointed by the tele- 
scope during the OFF-source runs as well as the date 
when each run was conducted. Each night was divided 
into two or three series of observation modes, such as 
"ON-OFF" , "OFF-ON-OFF" , or "OFF~ON" observa- 
tions. The ON-source run was scheduled to contain the 
meridian passage of the target, and the OFF-source run 
to follow the zenith angle distribution of the ON-source 
observation. On the average, the OFF source region had 
an offset in the R.A. of more than 30° away from the 
target. 

In the 2004 observation, the data of each telescope of 
T2, T3 and T4 were recorded with GPS time stamps, 
independently from the other telescope, when more 
than four photomultiplier (PMT) signals exceeded 7.6 
p. e. (photo-electrons) in each telescope. In the stage of 
offline analysis, the GPS time was used as a token to find 
coincidence events for the three telescopes. In the 2007 
observation, a t rigger circuit was emp loyed into the elec- 
tronics system (jNishiiima et al.ll2005() so that the events 
of more than two telescopes in coincidence were selected 
to be recorded. 

The total efficiency of light collection, including the 
reflectivity of the segmented mirrors, the light guides, 
and the quantum efficiency of the photomultiplier tubes, 
was monitored by "muon-ring analysis" of the Ch crenkov 
lights radiated by cosmic-ray muons (jEnomoto et al.l 
l2006a[ ). during the period of the observation. The de- 
terioration ratios of the reflectivity relative to the initial 
value when the mirrors were made were measured to be 
0.55, 0.60, and 0.75, for T2, T3, and T4, respectively, 
in 2004. The telescope T4 was built recently, and show 
the highest value of reflectivity. The deterioration rate 
was found to be slightly higher than the 10% level per 
year. Before the 2007 observation was commenced, mir- 
ror cleaning work using water was carried out. The re- 
flectivity was, however, not recovered to the initial value. 
The reflectivity of T2 became too low to be known as be- 
ing available from the calibration data of the muon ring. 
The two- fold coincidence data (i.e., T3 and T4) were used 
in the 2007 analysis. 

We selected those events in which the images of the 
Cherenkov light consisted of at least flve adjacent pix- 
els exceeding the 5 p.e. threshold (called the "cluster 
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TABLE 1 

Pointing positions of OFF source runs. "Offset" is the 

DECLINATION OFFSET FROM PSR B1706-44. THE DIRECTION OF 
THE CENTER OF FIELD OF THE VIEW OF THE TELESCOPE IS SHOWN 
IN THE EQUATORIAL COORDINATE OF R.A. AND 5, AS WELL AS IN 

Galactic latitude b and longitude I. 



Date 


Offset 


R.A. 


<5 


I 


b 


YYYY/MM/DD 


degree 


degree 


degree 


degree 


degree 


2004/Jul/ll 





243.89 


-44.48 


336.75 


4.57 


2004/Jul/13 





311.14 


-44.48 


356.17 


-38.35 


2004/Jul/14 





334.43 


-44.48 


353.19 


-54.79 


2004/Jul/17 





334.13 


-44.48 


353.29 


-54.58 


2004/Jul/18 





331.53 


-44.48 


354.10 


-52.79 


2004/Jul/19 





334.97 


-44.48 


353.00 


-55.16 


2007/Apr/14 


+0.5 


227.69 


-43.98 


327.79 


12.06 


2007/Apr/14 


-0.5 


227.69 


-44.98 


327.26 


11.20 


2007/ Apr/15 


-0.5 


213.00 


-44.98 


317.59 


15.58 


2007/ Apr/15 


+0.5 


213.00 


-43.98 


317.92 


16.53 


2007/ Apr/16 


-0.5 


212.30 


-44.98 


317.10 


15.73 


2007/ Apr/16 


+0.5 


212.30 


-43.98 


317.42 


16.68 


2007/ Apr/17 


-0.5 


206.97 


-44.98 


313.31 


16.74 


2007/ Apr/17 


+0.5 


206.97 


-43.98 


313.55 


17.72 


2007/ Apr/17 


+0.5 


204.97 


-43.98 


312.07 


18.02 


2007/ Apr/17 


-0.5 


204.97 


-44.98 


311.87 


17.04 


2007/ Apr/19 


-0.5 


209.80 


-44.98 


315.33 


16.25 


2007/ Apr/19 


+0.5 


209.80 


-43.98 


315.62 


17.21 


2007/ Apr/21 


+0.5 


197.26 


-43.98 


306.29 


18.77 


2007/ Apr/21 


-0.5 


197.26 


-44.98 


306.21 


17.77 


2007/Apr/22 


-0.5 


202.25 


-44.98 


309.89 


17.38 


2007/Apr/22 


+0.5 


202.25 


-43.98 


310.05 


18.36 


2007/May/12 


-0.5 


237.43 


-44.98 


332.96 


7.24 


2007/Mar/12 


+0.5 


237.43 


-43.98 


333.60 


8.02 


2007/May/15 


+0.5 


197.43 


-43.98 


306.41 


18.76 


2007/May/15 


-0.5 


197.43 


-44.98 


306.33 


17.76 


2007/May/16 


-0.5 


189.93 


-44.98 


300.76 


17.84 


2007/May/16 


+0.5 


189.93 


-43.98 


300.71 


18.83 


2007/May/17 


+0.5 


221.18 


-43.98 


323.56 


14.30 


2007/May/17 


-0.5 


221.18 


-44.98 


323.12 


13.40 


2007/May/20 


+0.5 


214.93 


-43.98 


319.28 


16.06 


2007/May/20 


-0.5 


214.93 


-44.98 


318.92 


15.13 


2007/May/20 


+0.5 


279.93 


-43.98 


351.48 


-16.50 


2007/May/20 


-0.5 


279.93 


-44.98 


350.51 


-16.87 


2007/May/21 


+0.5 


227.43 


-43.98 


327.62 


12.16 


2007/May/21 


-0.5 


227.43 


-44.98 


327.10 


11.30 


2007/May/21 


+0.5 


282.43 


-43.98 


352.16 


-18.18 


2007/May/21 


-0.5 


282.43 


-44.98 


351.18 


-18.53 


2007/Jun/13 


+0.5 


219.93 


-43.98 


322.72 


14.68 


2007/Jun/13 


-0.5 


219.93 


-44.98 


322.29 


13.78 


2007/Jun/13 


-0.5 


288.63 


-44.98 


352.64 


-22.69 


2007/Jun/13 


+0.5 


288.63 


-43.98 


353.68 


-22.41 


2007/Jun/14 


-0.5 


220.08 


-43.98 


322.83 


14.64 


2007/Jun/14 


-0.5 


220.08 


-44.98 


322.40 


13.73 


2007/Jun/14 


-0.5 


288.86 


-44.98 


352.69 


-22.85 


2007/Jun/14 


+0.5 


288.86 


-43.98 


353.73 


-22.57 


2007/Jun/15 


+0.5 


220.08 


-43.98 


322.83 


14.64 


2007/Jun/15 


-0.5 


220.08 


-44.98 


322.40 


13.73 


2007/Jun/15 


-0.5 


288.86 


-44.98 


352.69 


-22.85 


2007/Jun/15 


+0.5 


288.86 


-43.98 


353.73 


-22.57 



events"). For the 2004 runs, the frequency of the three- 
fold coincidence was S'^IO Hz, giving the rate of the clus- 
ter events as 6~7Hz. In 2007, the two-fold coincidence 
between T3 and T4 gave a trigger rate of 12^^14 Hz, 
which was reduced to 6~7 Hz for the cluster events. A 
cloud in the sky caused a low trigger rate when it blocked 
the direction of the telescope pointed, and events during 
the time of such a low-rate trigger were excluded from 
the analysis. The effective observation times for the 2004 
and 2007 observations were 996 (998) and 2187 (2386) 



minutes, respectively, where the time for an OFF-source 
run is shown in parentheses. The mean of the zenith an- 
gle of the data used for the analysis was 19.8° and 18.6°, 
respectively. 

4. METHOD OF ANALYSIS 

Since the ON-source run sees the Galactic plane and 
the OFF-source direction is away from the plane, as de- 
scribed in Table [U the amount of night sky background 
(NSB) hght differs between the ON- and OFF-source 
data. The region of the sky for the ON-source run is 
brighter than the OFF case. This may cause a spuri- 
ous effect of having excess events created by star light 
along the Milky Way, when we subtract the number of 
OFF-source events from ON-source ones. 

The NSB affects the distribution of ADC (Analog to 
Digital Converter) counts in the electronics circuit, and 
in turn, the degree of NSB is monitored in each obser- 
vation run. Darker NSB in an OFF-source observation 
can be "corrected for" to match the NSB in the case 
of ON-sou rce run, by ap plying the method of "software 
padding" (jCawlev 1119931 ) to add extra background light 
onto the photomultiplier tubes of OFF-source data. We 
applied a somewhat simplified software padding, when 
comp ared with the method developed by the Whipple 
group (fCawlev Ifl993h . Extra background lights of a com- 
mon shape of Gaussian fluctuation are added to ADC 
data of all PMTs for all OFF-source data in common. 
This padding procedure was applied only to the PMTs 
of having a signal of TDC (Time-to-Digital Converter) 
hit, since the CANGAROO-III analysis method requires 
TDC timing information. 

The effect that may be caused by adding more NSB as 
an artifact to the pixels of the photomultiplier tubes was 
studied by a Monte-Carlo simulation. It was not easy 
to precisely estimate and correct for the effects. How- 
ever, the spurious excess events could be made not to be 
effective by raising the threshold for discriminating the 
Cherenkov light signal in each pixel of the camera. We 
applied a threshold of as high as 8 or 10 p.e. for 2004 or 
2006 data, respectively, to each photomultiplier pixel of 
the clustering Cherenkov light image. 

The analysis procedures used in this work w ere basi- 
cally the same as those described in lKabuki et a l. (2007), 
but with m ore details to be foun d in lEnomoto et all 
([20G6a) and lEnomoto et all (|2006bD . As the fir st step, 
the Hillas parameters of a Cherenkov light image (jHillas I 
119851) were calculated for three or two telescopes. The 
arrival direction was determined by the condition of min- 
imizing the sum of the squared widths (xo! weighted by 
the photon yield) of the images seen from the assumed 
arrival position, in addition by putting a constraint on 
the Distance from the arrival position to the center of 
image. 

As a measure of gamma-ray likeliness of each event, we 
used the Fisher Discriminant (here after called as FD) 
(|Fisheilll936l : lEnomoto et al.l[2006ah . which is given by a 
linear combination of the input parameters, as 



The coefficients ai were chosen to realize the best 
separation between the gamma-ray and background 
events, where Pi is a component of vector, P ~ 
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{W2,W3,W4, L2, L3, L4) or P = {W3,W4, L3, L4), 
each of W2,W3,W4:, L2, L3, LA being the energy- 
corrected Widths or Lengths observed with tele- 
scopes T2, T3, and T4; for details see a reference 
(jEnomoto et al.l[2006a[ ) . We excluded from analysis those 
events in which the photomultiplier pixels in the out- 
ermost layer of the camera were hit by photons ("edge 
cut"), because the image of Cherenkov light of the events 
was deformed, particularly in the distribution of the 
Length parameter. 

We estimated the response function, f^, of FD for 
gamma-rays from a Monte-Carlo simulation, and ob- 
tained the function /{, for background events from the 
distribution of the OFF-source data. Note that the en- 
ergy spectrum proportional to E'' was used for gamma- 
rays in the Monte-Carlo simulation, where E is gamma- 
ray energy and 7=— 2.1. Then, the observed distribution, 
Eon, as a function of ED for ON-source data was equated 
to a linear combination of /-y and fb, as 



E„: 



Nbh + N^U 



Eh + E^ 



The number of background and gamma-ray events con- 
tained in the ON-source data were designated as Nf, and 
iV^, and were determined by fitting the ED distribution 
of the ON-source data to Ni,fb + N^f^ under the condi- 
tion of Non = Nb + N^, where Non is the total number of 
events of the ON-source data. The functions Eb — Nbfb 
and E^ = N^f^ represent the distribution of the back- 
ground and gamma-ray events with the total number of 
Nb and N^, respectively. 

The response function {f-y and fb) depends on its loca- 
tion in the field of view of the telescope, and thus a fitting 
for ON-source data by the sum of contributions from 
and fb was carried out for every arrival direction inves- 
tigated, within the field of view of the telescope. This 
way of inferring iV-y, the number of gamma-ray events, 
was adopted when we used the method (a). 

The method (a) of inferring by fitting the ED 
distribution was applied to the data of CANGAROO- 
III, to bring about the detection of g amma-ray signals 
successfully from point-like sources (jSakamoto et al.. I 
I2OO7I: iNakamori et al. II2007D as well as from a source of 
considerable extension (jEnomoto et al.ll2006bf l. and also 
to set upper limits on gamma-rays from o ther objects 
(|Kabuki et al.ll2007t lE^omoto et al."1l2007D . Note that 
with this method of fitting FD, it was not required to 
normalize the number of events by the observation time, 
in the case that the observation times of ON- and OFF- 
source runs were different from each other. 

A Monte Carlo simulation shows that gamma-ray 
events are likely to have ED larger than —0.5, which 
is consistent with the gamma-ray candidate events ex- 
tracted from the observation data by using method (a), 
as is presented in the next section. The events of ED 
larger than —0.5 were selected as gamma-ray enriched 
ones in the case of method (b): the wobble and ring 
background method. 

5. RESULTS WITH THE METHOD OF ON AND OFF RUNS: 
METHOD (A) 

The arrival directions in the equatorial coordinate 
(J2000) of observed events were sorted into 19 x 19=361 
cells with each having a 0.2° x 0.2° size in a field of 
view of ^FoV = 3.8° x 3.8°. For each cell of the arrival 



direction, the ED distribution of ON- and OFF-source 
runs was constructed to yield the parameter N^, i.e. the 
number of gamma-ray events, giving a map of gamma- 
ray-like events of the 2004 data, as plotted in Fig. [TJ 
Smoothing of taking the average of the neighboring 3 x 
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Fig. 1. — Map of excess events of 2004 data by using method (a) 
of the ON- and OFF source observations. The number of excess 
events per 0.2° x 0.2° cell is plotted in the equatorial coordinate. 
The number of excess events was smoothed by taking the average 
of the 3x3 = 9 cells around each investigated direction. The 
black cross at the center of the map indicated the position of the 
pulsar PSR B1706-44, the radius of the black circle represents the 
point spread function (PSF) of 59q = 0.24°, and the circle of white 
dotted line shows the region within 1° radius from the pulsar. The 
inserted white contours arc the hard-X-ray map obtained from the 
ROSAT satellite (iNASA n2007f) . 



3 cells was carried out for each investigated direction in 
the morphology map of the figure. The black cross at 
the center of the map indicates the position of the pul- 
sar PSR B1706— 44, the radius of the black circle shows 
the point spread function (PSF) of J^o — 0.24°, and the 
circle drawn with the white dotted line shows the re- 
gion within 1° distance from the pulsar. The white con- 
tours are the ROSAT h ard X-ray map taken from NASA 
Skyview (|NASA II2007I) . 

The excess, which is much larger than the statistical 
fiuctuation of ~ ±2 events per 0.2° x 0.2° area, shows 
a broader distribution than the point spread function. 
No prominent point-like excess appears either at the po- 
sition of the pulsar PSR B1706-44 or at any other di- 
rections in the vicinity of the pulsar, which is c onsistent 
with the H.E.S.S. observation (Aharonia n et al. ir2005a[ ) 
of no detection of point-like emission from the pulsar 
PSR B1706-44. 

The ED distribution of the observed events within a 
radius of 1.0° from PSR B1706-44 is shown in Fig. H 
The number of events of the ON-source data (the black 
cross) exceeds the distribution of the OFF-source events 
Eb (the green histogram) around ED ~ 0. The excess 
events, (Fo„— Ff,), indicated by the blue cross concentrate 
in ED = —0.5 to 0.5, is consistent with the gamma-ray 
distribution, Fy, inferred from the Monte-Carlo simula- 
tion and shown by the blue histogram. The right panel 
shows a magnified view of F^ and (Fo„ — Eb). The num- 
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Fig. 2.— FD distributions within 1.0° from PSR B1706-44, 
obtained from 2004 data by using method (a) of ON- and OFF- 
source observations. In the left panel, the data points with error 
bars (statistical errors) are the ON-source run, i.e., Fan by the black 
cross. The green histogram represents Fj,, the OFF-source data, 
the crosses of blue color are the excess counts of Fon — Fj,, and the 
histogram of red color shows Fj, the response of gamma-rays (fj) 
inferred from a Monte-Carlo simulation multiplied by N.y. The 
right panel shows a magnified view of the excess events. 




0.1 0.2 0.3 0.4 0.5 0.60.7 0.8 0.9 1 
(degree) 

Fig. 3. — Excess count map per solid angle (arbitrary unit) in 
6 space (the upper panel). Result of 2004 data by using method 
(a) of ON- and OFF-source observations. The middle panel is the 
estimated acceptance normalized at 6 at zero and the bottom is 
the acceptance corrected 6 distribution. The vertical axis is in 
arbitrary unit. 



ber of excess deduced is 504 ± 82 (6.1 a). 

The map of Fig. [1] might give an impression that the 
angular extent of the extended emission observed can be 
definitely determined to be as large as 1°. However, we 
have a difficulty to make a conclusive argument about 
the angular extent beyond 1°, since the "effective" field 
of view of our observation was limited to within about 
1° (half of the camera radius), due to the effect that 
a Cherenkov light image is deformed when a part of it 
is located near the outer edge of the camera. In order 
to estimate the angular extent of the excess events, we 
plotted in the top panel of Fig. [3] the distribution of 
gamma-ray-like events versus 9, where 9 is the angular 
distance of the arrival direction of each event from the 
black cross at the center of the map, the position of PSR 
B1706— 44. The detection efficiency as a function of the 
angle from the center of the field of view was estimated 
by a Monte-Carlo simulation, and also by using the num- 
ber of gamma-ray-like events of FD > —0.5 in an OFF- 
source run. As presented by the plot in the middle panel 
of Fig. [21 the efficiency decreases with 9, due mainly 
to the "edge cut". The acceptance of detection differs 
slightly between the 2004 and 2007 year observations, as 
shown by Fig. ID The acceptance for the 2007 data (the 
solid line) is wider than the 2004 one (the dashed line), 
due to the two-fold analysis for 2007 data when compared 
with the case of the three-fold analysis for 2004 data. 

The 9 distribution, after being corrected for the detec- 
tion efficiency, is shown in the bottom panel of Fig. [3l 
A correction for the acceptance was also applied to the 
morphology map of Fig. [U and is presented in Fig. [5l 
We limited the plot of data in this figure to be within 
a region, less than 1° from PSR B1706— 44. The accep- 
tance decreases with 9 rapidly outside of 1°, and the data 
seriously suffer from statistical fluctuation and system- 
atical error at 6* > 1° 




0.4 



0.2 - 



0.10.2 0.3 0.40.5 0.60.70.8 0.9 1 

(degree) 

Fig. 4. — Relative acceptance versus angular distance from the 
center of the field of view. The dashed line is for the three-fold 
analysis and the solid for the two-fold one. 

The result from data taken in 2007 is shown in Fig. 
[6l the FD distribution is in the left panel and the mor- 
phology map of the arrival direction is in the right, the 
number of gamma-ray-like events of the ON-source over 
the OFF-source runs is 627 ± 127 events for the region 
within 9 < 1.0° (5.1 a). 

The 6 plot from the combined data of years 2004 and 
2007 after the acceptance correction is shown in Fig. [71 
The x^/d.o.f value (the chi-square value per degree of 
freedom) calculated for flat distribution, i.e. the case 
of constant excess counts against 9, is 16.8/7. When 
we fit the distribution by a test function oc 1-9 /a from 
the pulsar position, the best fit for the observed dis- 
tribution was obtained when a — 1.36° ± 0.21° with 
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Fig. 5. — Map of excess events of 2004 data by using method 
(a) of ON- and OFF-source observations after an acceptance cor- 
rection. The definitions are the same with those in Fig. [T] 



ence in the ON/OFF distributions of the elevation an- 
gle. After a correction, however, this systematic error 
is reduced which was proved by the analysis of the off- 
Galactic-Plane data. The systematic error level is con- 
sidered to be less than 5%. The bright optical lights 
along the Galactic plane can deform the image parame- 
ters. We checked the FD responses by increasing NSB in 
the Monte-Carlo simulations and the OFF data (by soft- 
ware padding). When NSB increases, small sized images 
which look like gamma-ray images increase. These im- 
ages can be rejected by increasing the threshold of pixel. 
This is a reason why we used 10 p.e. threshold while the 
standard analysis used 5 p.e. We increased NSB amount 
and/or decrease this pixel threshold and check the differ- 
ences in the obtained gamma-ray fluxes. We, however, 
found that the main systematic error are due to NSB. 
The flat component observed in this analysis still remain 
to be statistically significant even in the worst case. The 
fluctuations in the observed flux changes approximately 
30% level. Our claim for the flat component is, there- 
fore, marginal. The resulting uncertainties in gamma-ray 
fluxes is -1-47, and -37%. 



TABLE 2 

Differential Flux from the area of 1.0° radius 
(corresponding to 9.57 X 10~'* sr) around PSR B1706— 44. 



Epoch <E> 
year TeV 



dF/dE 



10- 



^TeV- 



10" 



A(dF/dE) 

^cm^^g-iTcY- 



2004 



2007 



0.904 
1.00 
1.21 
1.05 
1.14 
1.29 
1.38 
1.54 



7.7 
5.3 
2.3 
3.1 
1.8 
2.1 
1.4 
1.5 



1.9 
1.0 
0.6 
1.1 
1.2 
0.7 
0.7 
0.3 



X^/d.o.f. = 3.7/6. In the case of fitting by a Gaus- 
sian function oc exp{—9^/{2a'^)), the x^/d.o.f value was 
4.0/6, with cr = 0.62° ±0.10°. Fitting was also made for 
the case of "constant value plus Gaussian function" , i. e. 
oc Nc+a-exp{-e'^/{2a'^)), yielding results ofiV^ = 19±6, 
a = 40 ± 11, and a = 0.34° ± 0.10° with xVd-O.f equal 
to 1.9/5. A morphology map for the combined data of 
years 2004 and 2007 after the acceptance correction is 
shown in Fig. [51 

The flux integrated within 1.0° from PSR B1706-44 
is plotted in Fig. [9] against gamma-ray energy E. The 
green points indicate the 2004 data and the black ones 
the 2007 data, respectively, with the error bars showing 
the statistical fluctuation. The flux of th e Cra b Nebula, 
reported by HEGRA (jAharonian et al. Il2000( l. is shown 
by the red line for a comparison. The flux is given nu- 
merically as a function of E in Table [2l 

In addition to our standard systematic errors such 
as on energy determination (photon collection efficien- 
cies), we investigated the possibility of systematic dif- 
ference of FD (gamma-ray likeliness) between two data 
obtained in different pointings. The width and length 
(image parameter) have a dependence on the elevation 
angle. The image becomes smaller in smaller elevation. 
Without any correction, the response of FD for back- 
ground events in the ON observation sometimes differ 
from that of the OFF. This is mainly due to the differ- 



6. RESULTS WITH WOBBLE AND THE RING 
BACKGROUND METHOD: METHOD (B) 

To the 2007 data, which were taken in the wobble 
mode observation, the "wobble" and the "ring back- 
ground method" can be directly applied in a way sim- 
ilar to what is used by the H.E.S.S. Group. The arrival 
directions along the circumference of its radius equal to 
9ri = 0.5° from the center of the field of view are guar- 
anteed to have the same, uniform detection efficiency. 
Along the circle of 9ri = 0.5°, six points that are suc- 
cessively at every 60° of the opening angle around the 
center of the field of view were selected as the central 
points of 6* = for calculating the 0^ distribution of Fig. 
[TUl where 6 is the angular distance from PSR B1706— 44. 
One of the six points corresponds to the center of "ON- 
source" directions, i.e. at the position of PSR B1706— 44. 
The other five points were utilized as the "off-source di- 
rection" for emulating background events. The events 
of having FD > — 0.5 were selected as gamma-rays 
enriched events, and then the numbers of events from 
"ON-" and "off-source" directions were compared for es- 
timating the flux of gamma-rays as a function of 9^ . We 
call this method of using the circle of radius 9ri = 0.5° 
as method (bl). 

Fig. [To] shows the 9^ distributions from the 2007 data 
by taking the events of FD > —0.5. When the cut 9^ < 
0.06 degree^ was applied to the observed 9^ distribution, 
as matched to the point-spread function, the excess was 
below 100 events, giving the upper bound of the flux 
roughly lower than the level of 20% Crab. 

To the 2004 data, method (bl) could not be di- 
rectly applied, because the points along the circumfer- 
ence of 9ri = 0.5° radius, on which the direction of 
PSR B1706— 44 is located, did not have a constant, uni- 
form efficiency of acceptance. It was necessary to correct 
for the varying acceptance by using a Monte-Carlo sim- 
ulation. The 9^ distribution, after the number of events 
was corrected for is shown in Fig. 1111 There appears 
to be no statistically significant evidence of point-like 
gamma-ray emission. The upper limit was set at the po- 
sition of the pulsar PSR B1706-44 to be less than 20% 
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Fig. 6. — FD distribution (left panel) and acceptance corrected excess count map (right panel) for the 2007 data by using method (a) 
of ON- and OFF-source observations. The notations are the same as in Fig. [2] and [5] 
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Fig. 7. — Acceptance-corrected 8 plot from the excess events of 
2004 and 2007 data combined. 
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Fig. 8. — Acceptance-corrected morphology map of excess events 
from 2004 and 2007 data combined. The definitions of contours and 
circle are the same as in Fig. [5] 

Crab flux. 

For constructing a morphology map, we used f 9 x f 9 
FD histograms of ON-source data, each of which was 
constructed from the events falling in a 0.2° x 0.2° area, 
within the field of view ofnpoV = 3.8° x 3.8°. The back- 
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Fig. 9. — Differential flux from the area within 1.0° radius (solid 
angle of 9.6 X 10"*) around PSR B1706-44. The green points 
are from the 2004 runs and the black from the 2007 runs. The 
errors are statistical ones. The re d line indicates the flux from 
Crab Nebula observed by HEGRA HAharonian et al. |[2000h . 

ground was estimated from those events, of which the 
arrival direction is at distance 0^2 from the investigated 
direction, namely within a ring-shaped region between 
9r2 = 0.40° and 0.55° from the investigated direction. 
We call this way of estimating background events method 
(b2) of the ring background. It is noted that this method 
differs slightly from "the H.E.S.S. method", in using dif- 
ferent values for some parameters such as the angular 
width. Background subtraction was then carried out af- 
ter an acceptance correction. A map of excess events is 
plotted in Fig. [121 The distribution of excess events is 
consistent with the statistical fluctuation, with no indi- 
cation of point-like source seen. The result of applying 
method (b2) to the 2004 data is presented in Fig. [121 
No prominent excess appears beyond the statistical fluc- 
tuation in the morphology map around PSR B1706-44. 

7. DISCUSSIONS 
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Fig. 10. — distributions of the 2007 data obtained by method 
(bl), the wobble method. The black crosses indicate the number of 
events obtained as the "ON-source counts" . The black histogram 
shows that obtained as the background number of events. The red 
points with error bars are the subtracted data of "ON-source" mi- 
nus the emulated background events, i.e. the "off-source" events. 
The green line indicates the zero position of the number of events. 
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Fig. 11. — 9 distributions of the 2004 data obtained by method 
(bl); the wobble method. The notations are the same as those in 
Fig. [To] 

The method (a) of having two separate observations 
of ON- and OFF-source runs has indicated emission of 
TeV gamma-rays from an extended region around the 
pulsar PSR B1706-44. On the other hand, method (b) of 
either the wobble or the ring background method, which 
is based on the subtraction of background events by using 
the neighboring region displaced by an angular distance 
of ~ 0.5°, could not confirm the result by method (a). 
The excess observed from method (b) is consistent with 
the statistical fluctuation, although arguments could be 
made about some spatial correlation that may exist in 
the morphological map of excess counts. 

The two differing results obtained through the two 
methods, both of which were applied to the same ON- 
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Fig. 12. — Morphology map of excess events for the 2007 data 
obtained by method (b2) of the ring background method. Details 
of the analysis are to be found in the text. 
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Fig. 13. — Morphology map of excess events obtained from the 
2004 data obtained by method (b2) of the ring background method. 
Details of the analysis are to be found in the text. 

source data, can, however, be understood as not in con- 
tradiction with each other, if we admit that gamma-ray 
events are contained in the dataset which is utilized for 
emulating background events. 

As can be seen in Fig. [21 the candidate events for 
gamma-rays are only ~ 2% of the total number of ON- 
source events. Such a small ratio of the signal-to-noise 
demands us to subtract from the ON-source data such 
events that are emulated from observed data containing 
presumably no gamma-ray events and are regarded as 
being equivalent to background events in the ON-source 
data. However, if a considerable amount of gamma-ray 
events are included in the dataset for emulating back- 
ground, the subtraction can diminish and even kill a 
gamma-ray signal that may exist. Such a chance may 
be more likely to occur in method (b) than in the case 
of method (a) of using the data from the region off the 
Galactic disk. 
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Diffuse gamma-rays from the Galactic disk should be 
subtracted to determine the genuine intensity of individ- 
ual sources, and thus the method (b) might be prefer- 
able to measure the flux from a source embedded in 
the Galactic disk. In fact, the diffuse gamma rays from 
the Galactic disk is expected to have a flux as large as 
Fdisk = 10"*^ ^ 10^^ cm-^s-iTeV^^sr^i at ITeV, when 
extrapolated from the 10 GeV energy region, or taking 
the c ase of the Galactic Center region (jAharonian et alH 
|2006() . Our result for the flux within an area of 1° ra- 
dius is 2.2 X 10~^^ cm~^s~^, as large as the intensity 
of the Crab nebula. The flux {Fdisk), when integrated 
over solid angle of a 0.2° x 0.2° = 3.8 x 10~^sr, amounts 
to 10"-^^ ~ 10~^^cm~^s~^TeV^^, which is equivalent 
to 0.01 ~ 0.1 of the flux from the Crab nebula. This 
intensity is similar to the flux from Galactic sources 
that the scan survey of the Galactic plane uncovered 
(jAharonian et al. Il2005b[) . 

If the emission from the Galactic disk is adequately 
uniform, or in other words, if the characteristic angu- 
lar scale of variance Odisk of the disk emission is larger 
than the angular width, 9foV, of the field of view of tele- 
scope and the angular size source of gamma-ray source, 
method (a) works to overestimate the flux by an amount 
of ~ i^(QsourceY' ' Fdisk- In this coutcxt, method (b) can 
be trusted to yield the "correct" flux, though it tends to 
be affected by neighboring sources to give a lower flux. 

When the wobble and ring background method (b) are 
utilized, the size of the gamma-ray emission, 9 source, is 
presumed to be less than the radius of the ring. The 
condition for method (b) being useful can be expressed as 
59q <IC 9source ^ 9^1 {oT 9r2) ~ 9foV- It might be said, in 
this context, that the angular resolution and the field of 
view of our observation is not very appropriate, namely, 
not good or wide enough, to clarify the morphological 
structure peculiar to the vicinity of PSR B1706— 44 by 
using method (b). 

In principle, and also in actual cases, it is neither 
easy nor clear to strictly distinguish the emissions of the 
Galactic disk and the individual objects from each other. 
The emission from the Galactic disk is likely to have 
a "granular structure" , which is caused by: the spatial 
distribution of energetic electrons with short lifetime by 
cooling, an irregular magnetic field of varying strength, 
complex distribution of the matter density of molecular 
clouds, the diffusion process of cosmic rays escaping from 
the acceleration site, and so on. Thus, the structure of 
the disk emission of gamma-rays is considered to have 
a variety of spatial scales, varying from place to place, 
depending on the local environmental conditions of the 
Galactic disk. 

The "granular structure" of the disk emission possi- 
bly lead us to erroneous conclusions, such as confusion 
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concerning the source position, incorrect estimation of 
angular size of the extended sources, or spurious detec- 
tion of point-like sources etc., which are influenced also 
by the performance of the telescope in use and, in par- 
ticular, by the scale size of the "granularity" 9disk- The 
energy spectrum of the obtained gamma-ray flux can also 
be affected by the gamma-rays in the background data. 

8. CONCLUSION 

Observations for about 50 hours with CANGAROO- 
III telescope system have given an indication of ex- 
tended emission of TeV gamma rays around the pul- 
sar PSR B1706-44. The strength of the signal de- 
pends on how we estimate angular size of the ex- 
tended emission. The total flux at 1 TeV is (4.7 ± 
0.7) X 10-"(£:/lTeV)-3 i±°-^ cm-^s-^TeV-i, when in- 
tegrated for incident angles within a circle of 1° radius. 
This corresponds to (4.9 ± 0.7) x IQ-^iE/lTeY)-^ '^^" '^ 
cm^^s^^TeV^^sr"^ in unit of "per solid angle". After 
integration of the gamma-ray energy, E, it is 2.2x 10^^^ 
cm~^s~^ for E >lTeV, which is as large as the Crab flux 
of 1.8 X 10"" cm-2s-i. 

The intensity within the area corresponding to the 
point spread function, 9 < 0.24° from PSR B1706-44, is 
(3.0 ± 0.6) X 10"^^ cm~^s~^ for gamma-ray energy E > 
1 TeV. The flux corresponds to 17% of the Crab flux at 1 
TeV, setting constraint on the emission from a compact 
source, which may underlie below the extended emission. 
The relative excess of this region compared with that of 
0.4° < 6* < 0.6° is 6 ± 4 % Crab. 

On the other hand, a statistically significant result 
could not be obtained from method (b) of the wobble 
and ring background analysis. The 2(T upper limit on 
the emission within 0.24° radius from PSR B1706— 44 is 
1.8 X 10"^^ cm~^s"^ at 1 TeV, which corresponds to 10% 
of the Crab flux at 1 TeV. 

The extended emission with method (a), together with 
the result using method (b), suggests complex struc- 
ture of TeV gamma-ray emission existing in the vicin- 
ity of PSR B1706— 44. A deeper investigation and fur- 
ther efforts for improving the technique of Imaging Air 
Cherenkov telescope remain to be pursued, in order to 
resolve and distinguish Galactic sources of TeV gamma- 
rays against the diffusive emission of the Galactic disk. 
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